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DIRECT COMPUTATION OF STRESSES IN LINEAR ELASTICITY 


WEIFENG QIU, MINGLEI WANG, AND JIAHAO ZHANG 

Abstract. We present a new finite element method based on the formulation introduced 
by Philippe G. Ciarlct and Patrick Ciarlet, Jr. in [Math. Models Methods Appl. Sci., 
15 (2005), pp. 259-571 ], which approximates strain tensor directly. We also show the 
convergence rate of strain tensor is optimal. This work is a non-trivial generalization of its 
two dimensional analogue in [Math. Models Methods Appl. Sci., 19 (2009), pp. 10)3-1064] 


1. INTRODUCTION 

This is a continuation of the two-part article proposed by Ciarlet & Ciarlet, Jr. [21 Hj. The 
main objective of this article is to introduce and analyze a direct finite element approximation 
of the minimization problem j(e) = inf j (e), which will be introduced in (12.51) . to compute 

e£E(fl) 

stresses inside an elastic body precisely in three dimensional space. The notations and ideas 
used in this paper to explain our finite element method are due to Ciarlet & Ciarlet, Jr. (See 

CD- 

Let S 3 denote the space of all 3x3 symmetric matrix. Let 12 be an open, bounded, 
connected subset of M 3 with Lipschitz boundary. Let a : b denote the inner product of two 
matrices a and b. Now consider a homogeneous, isotropic linearly elastic body with Lame’s 
constants A > 0 and /x > 0, with 12 as its reference configuration, and subjected to applied 
body forces, of density f G L s(12;M 3 ) in its interior and of density g G L3(12;M 3 ) on its 
boundary T. Given any matrix e = (e^) G 8 3 , Ae G § 3 is defined by 

Ae = A(tre)I 3 + 2ge. 


The classical way to solve pure traction problem of three-dimensional linearized elasticity is 
to find a displacement vector field u G R * 1 (12; M 3 ) which satisfies 

J(u ) = inf J(v), 


with 

J{v) 

For all v G H 1 ( 12; M 3 ), where 


- / A\/ S v: X/ S vdx — L{v). 
2 Vo 


(i.i) 


L(v ) = f ■ vdx + / g ■ vdT, 
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and 

Vsv = ^(V« T + V«) e S 3 ), 

denotes the linearized strain tensor field associated with any vector field v G R 3 ). 

If the applied body forces hold for the compatibility condition L{v) = 0 for all v G R(Q), 
where 

R(Q) = {v G H 1 ^; R 3 ) : y s w = 0 in 12} = (v = a + b A x : a, b G R 3 }. 

Then inf, ug #i(Q ; R3) J(v) > —oo. Besides, thanks to the Korn's inequality , the solutions of 
the minimization problem above exist (See [6j) and they are unique up to the addition of 
any vector field v G i?(12). 

Instead of finding vector field u directly, Ciarlet & Ciarlet, Jr (See [2]) put forward a 
new method to define a new unknown e which is a d x d symmetric matrix field in d- 
dimension (d = 2,3). It is proved by Ciarlet & Ciarlet,Jr in [2] that in both two and three 
dimensions (note that, Antman [Tj already proposed similar idea without proof in 1976 in 
three-dimensional nonlinear elasticity, while ours is concerned about linear elasticity), if 
e G E_{Q) where £7(12) = {e| curl curl e = 0; e G L s 2 { fl)}, then any vector field satisfying 
= e lies in the set {u|u = v + r for some r G _R(f2)} and the mapping n : E_(Vt) —)■ v G 
H 1 ( fl; R 3 ) = 7/ 1 (h2; R 3 )/R(Q) such that y s r = e is an isomorphism between AT 1 (C2; R 3 ) and 
E_(Q) (See HHIIJEIIE]). Since the relationship between e and v is clear (we choose e = S/ s v ), 
the minimization problem is converted to: 

]{e) = inf j(e) 

i(e) = ^ f Ae : edx-l(e), (1.2) 

for all e G E_(Q), where 

/(e) = L ok, 

and we can further show that £ = \/ s u. 

Let fibea triangulation of a polyhedral domain in R 3 and E_ h be a finite element subspace 
of We consider e h G E ,h to be 3 x 3 symmetric matrix with piecewise constant element. 

Our goal is to illustrate that how e h can satisfy the condition 

curl curl e h = 0 (1.3) 

Note that the similar ideas of the mixed finite element methods for linearized elasticity 
have been discussed in P [TUI jTTj, 121 031 03] • 

In order to achieve (1.3), the six degrees of freedom that define the elements e h G E h must 
be supported by edges (See Lemma 13.Ill and they must also satisfy specific compatibility 
conditions (See Theorem 13.2113.31 and 13.4p . Then the minimization problem can be accom¬ 
plished. The associated finite elements thus provide examples of edge finite elements in the 
sense of Nedelec (See |T5l 116] ). 

Finally, we consider discrete problem and it naturally comes to find a discrete matrix field 
e h G E h such that 



has unique solution (See Theorem IJ.lll . and we establish convergence of the method (See 
Theorem I4.2[) . 
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Note that this approach is in a sense the ’’matrix-analog” of the approximation of the 
Stokes problem by means of the divergence-free finite elements of Crouzeix & Raviart (al¬ 
though theirs are non-conforming, whereas ours are not). 

In the following part of the paper, the details of the approach and its validation will 
be discussed in section 2. Then the computational space is constructed in section 3. As it 
showed above, the final section will be the discussion of discrete problem and the convergence 
of the method. 


2. Three dimensional linearized elasticity 


The basic notations and conceptions are introduced in section two of [1]. Let Xi denote 
the coordinates of a point x G R 3 , let d\ := d/d Xi and dij := d 2 /dxidxj , i,j G {1,2,3}. 
Given a smooth enough vector held v, we define the 3x3 matrix held s/v := ( diVj ). In 
this paper, we hrst consider the domain which is open, bounded and connected subset of 
R 3 with Lipschitz-continuous boundary. Let fl be a domain in R 3 . Given any vector held 
v G R 3 ) as a displacement held, define 

V. v ■= ^(V^ T + V«) e L 2 ( Q) := L 2 (ft; S 3 ), (2.1) 

as its associated symmetrized gradient matrix held. Let 

R(Q) := {r €H 1 (n-,R 3 );S 7 8 r = 0} (2.2) 

denote the space of inhnitesimal rigid displacement helds. It is obvious that r G R(£l) if and 
only if there exists a constant vector c G R 3 and a matrix B satisfying B T = —B such that 


r = Bx + c, 

for all x = (x\,X2,Xz) G fl. 

It is discussed in [I] that the associated pure traction minimization problem 
J(u) = inf J(v), where J(v) = - 

veHHQ ;F Q ' Wo 


A Vs V : XJ s vdx - L(v). 


(2.3) 


with 


L(v) — / f ■ vdx + / g-vdT, 


has one and only one solution if it satishes the compatibility condition L(r) = 0 for all 
r G -R(O), which can be supported by Korn’s inequality. 

Instead of seeking the displacement held u in the classical approach, the intrinsic approach 
views the linearized strain tensor held e := \/ s u directly as the primary unknown. This 
changing of variables can be accomplished with the following results (See Theorem 2.1,2.2 
and 2.3), which were proved in [2] by Ciarlet and Ciarlet,Jr. 


Theorem 2.1. Let Q be a simply-connected domain in R 3 and let e = (e^) G L 2 (fl) be a 
tensor field that satisfies 

curl curl e = 0 in R 3x3 ). 

Then there exists a vector field v G 7L 1 (h2; R 3 ) such that \7 s v = e in L 2 (Ll), and all the other 
solutions v are of the form 

v = v + r for some r G R(Q), 
where R(fil) is the space defined in (12.211 . 
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Theorem 2.2. Let D be a simply-connected domain in R 3 . Define the space 

E_(Ll) := {e G Z^(Q); curl curl e = 0 in JW 2 (12; R 3x3 )}, (2.4) 

and given any e G let v = n(e) denote the unique element in the space H l (Vt\ R 3 ) 

that satisfies e = \/ s v (See Theorem 2.1). Then the linear mapping 

k : E(Vt) H l ( Q;R 3 ) 

is an isomorphism between the Hilbert space E_(Ll) and R 3 ). 


Theorem 2.3. Let Ft be a simply-connected domain in R 3 . Then the minimization problem: 
Find e G _E(12) such that 

j(e) = inf j(e), where j(e) — ~ [ Ae:edx — l(e), (2.5) 

ee£(D) ' 2 J n 


with 


has one and only one solution e and 


/(e) = L on, 


£ := 

where u is the unique solution to problem (I2.3[i . 

3. Finite element space for strain tensor 

Same as Ciarlet & Ciarlet, Jr. did in [T] , we first describe a triangular finite element, 
which provides edge finite element. The length element is denote dl. 

Lemma 3.1. Let T be a non-degenerate tetrahedron with edges Si, 1 < i < 6. Given any edge 
Si ofT, let r* denote a unit vector parallel to Si, and let the degrees of freedom di, 1 < i < 6, 
be defined as 

dfie) := ! F-efidl VeG? 0 (T;S 3 ). (3.1) 

J Si 

Then, the set {dp, 1 < i < 6} is Po(T-,E> 3 )-unisolvent, i.e. a tensor field e G Po(T-,E> 3 ) is 
uniquely determined by the six numbers dfie), 1 < i < 6. 

Proof. The proof will be the same as that of Theorem 3.1 in [1]. □ 

From now on, denotes a polynomial domain in R 3 , and we consider triangulations T h 
of the set 0 by tetrahedrons T GT fe subjected to the usual conditions; in particular, all the 
tetrahedrons T G 7 h are non-degenerate. 

Given such a triangulation 7 h of 0, let Tfi denote the set of all “interior” edges found in 
7 h (i.e. that are not contained in the boundary <9f2), let Eg denote the set of all “boundary” 
edges found in 7 h (I.e. that are contained in <9f2), let A h denote the set of all “interior” 
vertices in 7 h (i.e. that are contained in 12), and let Aq denote the set of all “boundary” 
vertices in 7 h (i.e. that are contained in <9f2). 

We also assume that each interior or boundary edge cr G E^U Eg is oriented. For instance, 
if {a- 7 ; 1 < j < J} denotes the set of all the vertices found in 7 h and a = [a 1 , a- 7 ] G E ft U Eg 
with i < j, one may let r : |a 7 — a*| _1 (a 7 — a 1 ). 
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Theorem 3.2. Given any triangulation 7 h of Cl, define the finite element space 

if := {e h e L*(fi); e h \ T G PfiT- § 3 ) for all T G 7 h , and for all a G £\ (3.2) 

r ■ ( e h \ T )rdl is single-valued for all T G 7 h with a C T}. 

Then each tensor field e h G if is uniquely defined by the numbers d a (e h ),a G T, h UT,g where 
the degrees of freedom d a : E —* R are defined by 

d a (e h ) —It- (e h \ t)tcII for any T G 7 h with a C T. (3.3) 

J a 

Define the finite element space 

E h :={VsV h eh 2 s (nfiv h e V h }, (3.4) 

where, the space R(Cl) := {v = a + h A x; a, b G M 3 }, 

V h : = V h /R(Cl) with V h := {v h G C°(C2); v h \ T G PfiT- M 3 )}. (3.5) 

Then, 

c E h . (3.6) 

Besides, 

dimE = dim V h — 6. (3.7) 

Proof. According to Lemma ITT) each tensor field e h G E is uniquely defined by the numbers 
d a (e h ) in (13.3p . The proof of (13.6113.7p will be similar to that of Lemma 3.1 in [I], □ 

Lemma 3.3. Given any triangulation 7 h of Q, for any vertex a found in 7 h , let {T;T G 
^(a)} denote the set formed by all the tetrahedrons of7 h that have the vertex a in common, 
and let f2 a := Utst' 1 T . 

For any given vertex a found in 7 h , We denote N number of vertices in f2 a; A number of 
edges in Cl a . If a G Tfi, then 

A = (3N — 6) + (N b — 3), (3.8) 

where, Nj, is number of vertices in the boundary ofCl a . If a G Tig, then 

A = (3N - 6) + N ib , (3.9) 

where, Ni b is number of vertices in dCl a \dCl. 

Proof. If a G Eg, we denote A ib number of edges found in dCl a \dD. Then, 

A = A ib + 2 (N — N ib — 1) + N ib . 

And, we utilize (3.12) in [I] for dCl a \dCl, then 

A-ib + (A" — N ib — 1) = N ib + 2(.V — 1) — 3. 

So, we have that 



A — A ib + 2 (A" — Ni b — 1) + Ni b 
= [N ib + 2 (N — 1) — (N — - 2] + 2 (N - N ib - 1) + N ib 

= (3 N — 6) + N ib . 
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This proves (13.9p . 

If a £ E h , we take T £ 7 h (a) arbitrarily. Then, it is easy to see that A and N are the 
same number of edges and number of vertices in Q a \T, respectively. We treat dfl a \dT the 
same as dfl a \dfl in above argument. This will immediately prove (13. 8 j) . □ 

Theorem 3.4. Given any vertex a £ Ag U A h , there exist m-many linearly independent 
linear forms (p a j : E —> E (1 < j < m) such that 

curl curl e h = 0 in D'(12) for all e h £ E h , 
where the space E h is defined by 

E /l := {e h £ ff; < p a ,j(e h ) = 0,1 < j < m for all a £ Ag U A h }. (3.10) 

Here, m = N lb if a £ Ag, and m = N b — 3 if a £ A h . N ib and N b are defined in Lemma UTS 1 
More specifically, the coefficients of each linear form <p aj are explicitly computable func¬ 
tions of coordinates of the vertex a and the vertices of the tetrahedrons of 7 h that have a as 
a vertex. 

In addition, if the domain 12 is simply-connected, then E h = E . 

Proof. Given any vertex a £ Ag U A h , we denote Cl a the closure of union of all tetrahedrons 
which have vertex a. It is easy to see that 

dimE' |q o = number of edges in Cl a , dimE |n a = 3 x (number of vertices in Q a ) — 6. 
Then, by Lemma [3.31 

dimE^lf^ — dimE |^ o = N ib if a G Ag, dim e !'\— dimE |fj a = N b — 3 if a G A h . 

This implies that there exist m- many linearly independent linear form (p a j : E |q o —>■ ffi. 
(1 < j < m) such that E |n a is equal to 

{§. h e Pa,j(e h ) = 0,1 <j< m}. (3.11) 

Now, we take e h e E h arbitrarily. If (p a g(e h \n a ) = 0 for all 1 < j < m, then we have that 
er\ cia ^ K I ha- This implies that 

curl curl e h = 0 in D , (int(r2 a )). 

We notice that hi = UaeAgu A h i n f(^a)- So, we can conclude that 

curl curl e h = 0 in D'(12). 

And, for any e h £ E ? * and any a £ Ag U A h , 1 < j < m, we can define 

Pa,j(§_ ) ■— ^ Pa,ji,§L |n o )- 

So, we can conclude that 

curl curl e h = 0 in 2)'(H) for all e h £ E h . 

Obviously, if C E h . 

If the domain 12 is simply-connected, then using the same argument in the last paragraph 
of the proof of Theorem 3.3, we an conclude that E h = E . □ 
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4. The discrete problem and convergence 

In the following, Let II be a simply-connected polygonal domain in M 3 . The discrete 
problem is defined as the minimization problem (See (4.1)) below. 

Theorem 4.1. Given any triangulation I\ of PL, let E h be the finite element space defined 
in \3.1C\) . Then there exists one and only one £ h G ~E h such that 



where j is defined in (12.51) . Besides, let J be the function defined in (12.3f) . then e h = \/ s u h , 


where ii h is the unique solution to the minimization problem 

J{ii h ) = inf J(v h ). (4.2) 

■v h GV h 

Proof. The proof is the same as Theorem 4.1 in jl]. □ 

Finally, we discuss the convergence of the method. 

Theorem 4.2. Consider a regular family of triangulation of Pi. Then 

||s — —y 0 as h —y 0. (4-3) 

If u G H 2 (Pl\ M 3 ), there exists a constant C independent of h such that 

||li — — C\\u\\ H 2(n)h. (4.4) 

Proof. The proof is the same as Theorem 4.2 in jl]. □ 
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